Quantum corrections and bound-state effects in the energy relaxation of hot dense 

Hydrogen. 
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Simple analytic formulae for energy relaxation (ER) in electron-ion systems, with quantum cor- 
rections, ion dynamics and RPA-type screening are presented. ER in the presence of bound electrons 
is examined in view of of recent simulations for ER in hydrogen in the range 10^''-10^'' electrons/cc. 
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Introduction - The problem of energy relaxation (ER) 
in hot clcctron-ion systems arises in astrophysics, fusion 
physics, condensed-matter physics and chemistry [l|, |^. 
The large difference in masses Trig, mi between elec- 
trons and ions simplifies the problem somewhat since 
quasi-equilibrium systems, with an electron temperature 
Tg, and an ion temperature Ti can occur. Such two- 
temperature plasmas exist in material systems as diverse 
as warm dense matter 0], ultra-cold plasmas [J|, hot 
semiconductors |5| , and dense deuterium [&\ . The strong 
temperature dependence of thermonuclear processes im- 
ply that estimated burn rates depend sensitively on the 
accuracy of the temperature relaxation (TR) theory used. 
Thus a number of recent studies [l|, 0, Q have examined 
various aspects of two-temperature dense hydrogen and 
ER in regimes of densities and temperatures relevant to 
inertial confinement fusion (ICF)[3]. 

The earliest theories of ER in plasmas are due to Lan- 
dau [l3l and Spitzer [ll| (denoted L-S). The L-S ap- 
proach is applicable to weakly interacting fully-ionized 
plasmas in the classical regime. It is the Rutherford 
Coulomb scattering formula applied to Maxwellian distri- 
butions of ions and electrons. Analyses using the Fermi 
golden-rule (FGR) and coupled-mode (CM) extensions 
were given by Dharma-wardana et al. [l|, Hi Ej KL4| . 
Fokker-Plank approaches [13], kinetic-equation methods 
based on the relaxation of the one-particle distributions 
[isl [Tq . methods based on expansions in the couphng 
constant I17| or related techniques [l^, have been explored 



in recent studies [TT]. Molecular dynamics (MD) was used 
by Hansen and McDonald, (denoted HM) |l0| , using clas- 
sical potentials which incorporate the cutoffs used in the 
Coulomb logarithm of L-S theory. The recent studies of 
Refs. [l|, Q report more careful applications of the HM 
approach. The FGR and CM approaches have been ex- 
perimentally tested in semiconductor plasmas, but no ex- 
perimental results of TR are as yet available for hot dense 
plasmas. 

In fully-ionized hydrogen plasma the particle charges 
Zp, Ze are ±1. The mean electron- and proton densities 
n and p are identical. The ratio of a typical Coulomb 
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energy to the kinetic energy becomes, in the classical 
regime, F — \/{rgT), where T is the temperature in en- 
ergy units, and r^ — [3/(4?™)] ' is the radius of the 
Wigner-Seitz sphere of an electron or a proton. The 
properties of partially degenerate plasmas require two 
independent parameters, e.g., both F and 9 — T/Ep, 

where Ep — (Svr^n) /2 is the Fermi energy. Thus the 
regime of densities and temperatures studied in Ref. [8| 
involve, at one extreme, r^ = 25.25, 9 — 127 at Te=10 
eV., for rie — 10^" electrons/cc, while the another ex- 
treme is rs = 1.172, and 9 = 0.274, where Ue = lO^'' 
electrons/cc at T^—IQ eV. The latter is a significantly de- 
generate plasma where the validity of classical-simulation 
methods is suspect. The system at r^ ~ 25, rie = 10^° 
electrons/cc at T^ 10 eV contains Is to 3d bound states, 
but the system is essentially ionized since the bound- 
state occupations are negligible. On the other hand, 
the plasma at Ue — 10^^ electrons/cc, T^ = 10 eV, i.e., 
rs=5.441 is ~ 82% ionized at T^ = 10 eV and carries a Is 
bound state at an energy of -0.458 a.u., when calculated 
using the atom-plasma codes implemented by Perrot et 
al[20'|. The effect of such bound states cannot be included 
in the classical simulation method of Ref. 8] . 

In this communication we use the FGR approach and 
derive a simple analytic formula inclusive of leading 
quantum corrections, screening and ion-dynamics for for 
TR in the fully ionized limit, and compare it with the re- 
sults from recent classical simulations [8|. We then con- 
sider the effect of bound states, since they can have a 
significant effect on energy relaxation. Many complex 
processes become possible, but the system with a single 
hydrogenic boundstate is a useful case for sharpening our 
understanding of this relatively unknown regime where 
no previous results are available. 

Quantum transition rates.- Assuming that Tg > Ti to 
be specific, ER occurs via energy transfer from the ex- 
cited modes of the electron sub-system to the cold modes 
in the ion subsystem. The spectrum of the modes of the 
species j is given by the spectral function Aj{q,uj,Tj). 
These spectral functions are given by the imaginary 
parts of the corresponding dynamic response functions 
X'{k,uj), e.g, Eq. (16) of Ref. [H. The ER rate evalu- 
ated within the Fermi golden rule, Rfgr can be expressed 
in terms of the response functions of the plasma as given 



in Eqs. (4)-(7) of Ref.[3, and Eq. (15) of Ref. [Tj 
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In the above 6E/6t is the rate of change of the energy 
of the system, for time steps St significantly greater than 
the equihbriations times Tp, Tp which cstabhsh Tg and Tp 
of each subsystem. The relaxation of the whole system is 
determined by r^p such that Tep >> Tp > Te- For brevity 
we write 6E/St as dE/dt. we have used the spherical 
symmetry of the plasma to write scalars q, k instead of 
(f, q to simplify the notation. The non-interacting re- 
sponse function xo{1i^tT) at arbitrary degeneracies was 
given by Khanna and Glyde[2l|, and are used here in 
the generalized RPA x^(g,a;,T) form, j — e,p where, 
for example, temperature-dependent local field correc- 
tions Gee{k) may be included z2]. If Tg, Ti are both suffi- 
ciently large that AB — > 2{Te — Ti)/u!, and if the electron 
chemical potential fie < 0, useful analytical approxima- 
tions become available. The possibility of extracting a 
temperature-relaxation time Tgp from the relaxation rate 
exists only in this regime. Neglecting interactions, E 
becomes the kinetic energy. Using non-interacting clas- 
sical forms for 3xo(fcj'^) in Eq. [T] we obtain the well 
known Landau-Spitzer (L-S) form for the temperature 
relaxation time r, viz., 
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FIG. 1: The Landau-Spitzer Coulomb Logarithm (£) is com- 
pared with the classical MD simulations p] and the analytic 
Coulomb factor Q, Eg. 1161 derived from the FGR result. 
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Here cop^ is the plasma frequency of the species j = e,p. 
The effective temperature and the effective mass of the 
colliding pair are T^p and mep, with Tj in energy units. 
C is the "Coulomb logarithm". It depends on kmin 
and kmax, i-e., momentum cutoffs (or impact parame- 
ters) used for modeling the unscreened Coulomb collision. 
If interacting response functions (e.g, RPA) are used, 
single-particle modes become replaced almost completely 
by plasmon modes, and the interactions become dynam- 
ically screened. Any type of "static screening" must sat- 
isfy the /-sum rule if false results are to be avoided. The 
safest procedure is to do the w-integration imposing the 
/-sum rulejl4|]. Then ion dynamics are automatically 
preserved, and Eq. [1] simplifies to: 
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The anti-symmetry of the imaginary part has been used 
in the above. The electron-screening wavevector kgc at 
any degeneracy is obtained from the small-g limit of the 
finite-T Lindhard function. The k -^ 0-local field correc- 
tion, Gg^ at arbitrary degeneracy [22| can also be included 
in ksc via the following definitions. 
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However, in this study we find it sufficient to use he 
Debye-Hiikel form ksc = k'^n l^r electron screening. In 
approximating '^dx^'^ / doj we retain terms up to second 
order in Ti as displayed explicitly below: 






(— ) 



3/2 



2 2ria; ^r ■;£_££*£ j sinh(?iaj/2rt 



Trfc 



Tiuo/2Te 



Then Eq. [7] can be reduced to the form: 
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The exponential integral [23|, Ei{x) of Eq. [TH is eval- 
uated numerically via standard subroutines. Thus we 
see that the "Coulomb factor" Q is exactly analogous to 
the "Coulomb logarithm" of Eq. [6l but without ad hoc 
cutoffs. Q contains leading-order quantum corrections, 
ion-dynamics and electron screening. The expression for 
Q should be compared with a similar expression given by 
Brown et alpTf. 



^bps = 2 [log(l/Pe) 



■7 



(17) 



where 7=0.5772 is the Euler constant. At high Tg, this 
result approaches the L-S form more rapidly than Q. 

The interactions between the ion- and electron modes 
lead to ion-acoustic modes (coupled modes) [1]. The H- 
plasmas treated here may be considered relatively weekly 
coupled plasmas, and the correction from coupled- mode 
effects will be neglected. The results shown in Fig. [T] 
suggest that the classical potentials etc., used for Ts — 
5.44 and Ts — 25.25 need reconsideration, since these 
are weakly coupled plasmas where FGR methods should 
give good agreement. As seen from Fig. [21 the numerical 
results from Q and Cbps very similar. The case r^ — 1 
has been discussed in more detail elsewhere (l|. 

Energy Relaxation in the presence of bound states- The 
nuclei in the H-plasma with 10^^ electrons/cc, r^ — 5.44 
carry a single Is bound state of energy ~0.46 a.u. with 
the effective charge Z varying from 0.82 at Te=10 eV, to 
~ 1.0 at Te > 200 eV. If we use an "average ion" picture, 
we have a gas of ions with change Z, a boundstate occu- 
pation nis — 1 — Z electrons per ion, and the remain- 
ing ("free") electrons distributed in continuum states. 
The electrons in the bound states have an effective mass 
niep ~ rUe = 1, and hence equilibrate rapidly with the 
hot electrons in the continuum. That is, the distribution 
nis, and hence the degree of ionization Z is determined 
by the hot electron temperature Tg- (See fig.[2]for a plot 
of 1 — ^). The center of mass (CM) motion of the ions 
and bound electrons is determined by the kinetic energy 
of the mass Mqm = Mp + rUe — Mp, i.e., CM motion is 
at the temperature Tp. The energy exchanging collisions 
are between the hot electrons and the ions. The bound 
electrons (in their Is states) do not interact with their 
own binding nuclei. The energy mismatch in the spectral 
functions ensures that they do not significantly interact 
with other ions carrying bound states. This "average- 
ion" picture holds if the temperature relaxation time r is 
sufficiently long compared to the ionization-equlibriation 
time Tz- This picture is applicable to the system dis- 
cussed in this study, and the only source of energy relax- 
ation is via Coulomb collisions. 

An alternative picture holds if r is comparable to tz- 
Then the plasma consists of nonionized atomic hydrogen 
with fully occupied Is states at some temperature Ta of 
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FIG. 2: The effect of bound states in the H-plasma with 
Vs = 5.44, shown as Qb, compared with calculations without 
bound-state effects. Lbps is also shown to closely approxi- 
mate Q. The bound fraction 1 — Z, where Z is the degree of 
ionization, is shown in the inset. 



the atoms, free electrons at Te and protons at Tp. The 
atom temperature Ta has to be self consistently deter- 
mined via a theory of ionizing collisions as well as the 
Coulomb collisions discussed so far. This regime occurs 
for low Te,Tp and is not considered here. 

The energy E^ of the electron subsystem, and the 
energy Ei of the CM-ion subsystem in the average-ion 
model with hydrogenic states e„, occupations n„ (i' — 
n, Z, m being the set of quantum numbers of the state) is 
given by: 

Ej = 3nTp/2 + i?,L (18) 

Ee = 3nZTe/2 + n(l-Z)^(e^-/Xa)n^ (19) 



3nZTe/2 + n{l - Z)\is + E^.^ 



(20) 



The terms Ef^f. involve contributions from particle in- 
teractions, embedding effects, continuum lowering effects 
etc., which have been discussed at length in Ref. [24| . 
These effects can be neglected for the weekly-coupled H- 
plasmas studied here. The last equation limits the prob- 
lem to one bound state. The effects arising from the 
chemical potential /i^ = ^e{Te) + fJ^piTp) + Spa will be 
neglected in all the species, free and bound, as they are 



found to be small. Then the time dependence dEj/dt can 
be replace by dEj/dTj){dTj/dt) and we attempt to con- 
struct a formula for temperature relaxation. Unlike in 
simple classical plasmas, it is no longer possible to define 
a temperature relaxation time r. Hence, to be specific 
we may thermostat Ti and consider the relaxation rate 
dTe/dt, with {Z ^ 1) and without {Z — 1) bound-state 
effects. Thus we have: 
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This result shows that the effective ionization Z enters 
approximately as Z^ in the prefactor. The dependence 
of the ionization temperature, viz., dZ/dT^ was obtained 
from a DFT calculation[20]. ft may be argued that some 
of these corrections are more part of the specific heat 
of the electron system inclusive of bound states, and do 



not belong to the collision dynamics. However, merely to 
compare the Z = I, Z ^ I systems, we define 



Qb = {Z^/Bpe}Q{Zpe 



(24) 



We display Q and Qf, using the Eqs.[24l in Fig.[2l together 
with the classical MD, Chps and L-S results which treat 
only the fully ionized state. The inset to the figure shows 
the degree of association [1 — Z) as a function of ternper- 
ature, calculated using density- functional methods ,20|. 

Conclusion. -^e presented analytic formulae with 
quantum corrections, ion dynamics and screening . They 
have been applied to weakly coupled H-plasmas (rg = 
5.44, 25.25) where the physical approximations are felt 
to be reliable. The differences between the classical MD 
simulations and our results suggest the need for further 
review of the classical simulations as well as the physi- 
cal theory. The effect of the bounds states on a quan- 
tity nominally similar to the Coulomb logarithm has also 
been presented for the H-plasma at the density \Q^^ e/cc. 
This too is an area which requires further study and poses 
a novel challenge to simulations.. 
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